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Abstract 

We present an approach to derive a relativistic kinetic equation of the Vlasov type. Our 
approach is especially reliable for the description of quantum field systems with many internal 
degrees of freedom. The method is based on the Heisenberg picture and leads to a kinetic 
equation which fulfills the conservation laws. We apply the approach to the standard Walecka 
Lagrangian and an effective chiral Lagrangian. 
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1. Introduction 

Experiments at CERN and BNL are designed to investigate, the formation and decay 
of a quark gluon plasma in ultra-relativistic heavy-ion collisionsa. The very early stage of 
such a collision is dominated by non-equilibrium effects. The expected quark gluon plasma 
undergoes a transition to strong interacting and highly correlated hadronic matter if the 
system cools down and becomes dilute. The description of such a dynamical evolution of a 
very complex, relativistic system with many internal degrees of freedom is one of the most 
challenging problems in modern particle and non-equilibrium physics. 

Relativistic kinetic theory and hydrodynamics have been successfully emplov.ad for the 
description of non-equilibrium states of matter at high density and temperaturddtlu. These 
methods have been most-c^CiOijjiiiely applied in relativistic nuclear physics both in the in- 
termediate energy regioiuauOLTLJ, where the quark-gluon degrees of freedom arie-istill not 
relevant, and in the high energy domain of a non-equilibrium quark-gluon plasmaliil. In the 
construction of the relativistic kinetic theory, different approaches have been used. Within 
the contour Green function technique, relativistic transport equatieps of the Vlasev type have 
been developed by Kadanoff and BayrnES, Martin and SchwineepJ and Keldyshufl. That ap- 
proach has been studied formally by a large number of authorsu; and others have explored it 
in quantitative model studies. 

However, the theory can be developed further. The derivation of a relativistic transport 
equation that holds for systems with many inter»a|l degrees of freedom as well as the inclusion 
of collision and particle production source termscS are of particular interest. 

We investigate an approach based on the non-equilibrium statistical operatepiji Heisenberg 
dynamics that can be considered as a generalization of the Zubarev methodcj'tij to the rela- 
tivistic domain. We restrict ourselves to the mean field approximation but take into account 
the internal degrees of freedom of the constituents. 

For demonstration we use,two different models. The Walecka model as the standard version 
of quantum hadrodynamicstJ serves for the description of hadronic degrees of freedom of rela- 
tivistic nuclear matter for laboratory energies less than a few GeV/A. A further application is a 
chiral Lagrangian that is of particular interest because the concept of chiral symmetry breaking 
at large temperatures and density is supposed to be the driving force for the hadronisation of 
matter. Therefore microscopic quark models, such as the Nambu-Jona-Lasinio (NJL) model, 
may be useful for the dynamical description of the expected quark-hadron phase transitions. 

The paper is organized as follows. In Section 2 we present an overview of the method. 
In Sections 3 and 4 we demonstrate the derivation of the relativistic kinetic equation of the 
Vlasov type using the Walecka model and the NJL model, respectively. We summarize the 
results in Section 5. 

2. Equation of motion for the Wigner function 

In quantum field theory the equations of motion within the Heisenberg picture are given 

as 

id»A(x) = \A{x) 1 P»} , (1) 
where A(x) is an arbitrary local operator and is the total 4-momentum of the system, 

P" = J da,{n)T"" . (2) 

Here da M (n) denotes a vector element of an arbitrary space- like hyperplane with a time- like 
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normal vector n M (n 2 = 1), T ,J, "(x) is the energy-momentum tensor. For the fermionic sector 
the energy-momentum tensor T M " reads 

T»» = <?<</> + fyf d" - gTC{x). (3) 

The interaction is not specified yet and will be considered in the Walecka model (Section 3) 
and the NJL model (Section 4). 

Let us now consider the situation when the time-like direction in Minkowsky space is 
determined by an external condition. We suppose in the following that this direction is defined 
by the unit vector n M which fixes simultaneously the orientation of the hyperplane c(n) cil 
in the formulae (Q) and (|fl|). For the description of the dynamics of the system along the 
time-like direction n M and along the independent space-like directions on the hyperplane <j(n) 
, the boost transformation with "velocity" can be used for the equations of motion (hi). 
Technically, it is achieved with the help of the projection of Eqs. (|l|) on the direction n M and 
on the hyperplane cr(n). Convolution of (|l]) with n M leads to the dynamical equation 

i^§^ = [A(x),H(r)]. (4) 

The parameter 

r = n M x M (5) 

plays the role of a proper time in a new coordinate system. The derivative along the direction 
of n M in Eq.(0) is given as 

9 9 («\ 

a7 = n ^ ' (6) 



H(t) is the Hamiltonian of the system I 

H(t) = n M P" = / da{n)n^ v n v . (7) 



The remainder equations 

i-^A(x) = [A(x),H^r)] (8) 

intend for a description of infinitesimal transfers of a system on the space-like hyperplane a(n). 
The space-like derivative is defined by the relation 

— =K— , (9) 

dxj: v dx v ' 

where A M " = g>*" - n^n v is the projection operator, A^n,, = n M A Ml ' = 0. And finally, 

IP(t) = KP V = / da{n)n x T Xv K (10) 



is the space-like momentum vector of the system. H(r) and n M (r) are the generators of local 
infinitesimal normal and tangential diffeomorphismes to a hyperplane of a constant proper 
time. Note that a similar covariant decomposition of the gewen-ator of a motion group can be 
found by the quantization of the Einstein gravitation theoryca. 
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In order to derive the drift integral of the relativistic kinetic equation it is convenient to 
introduce the anti-commutation relations for fermionic fields: 

{1>(x),$(x')} = -iS(x-x'), Wi)^(/)} = {^(i),^')}=0 (11) 

for any points x* and x'* which belong to an arbitrary space-like hyperplane a . Let us note 
the following property of the commutation function of spinor fieldsB 

Si?) Uo =0 =n 5 (3) W , (12) 

and 

da(n\x)S a f3(x — x') y(x') = in^api/ix), x£a{n), (13) 



where y(x) is an arbitrary function of field operators and the spin is denoted by the indices a 
and /3. 

The standard definition of the one-particle covariant Wigner function of the Fermi subsys- 
tem reads 

f«0(x,p) = J dye-' py < P a p(x,y) >, (14) 

where < ... >= TrL.. ^} denotes the operation of statistical averaging with the single particle 
density matrix p HZrES in mean field approximation within the Heisenberg representation with 

P a p(x,y) = $p{x + y/2)i, a {x - y/2) . (15) 

Performing the r - differentiation of the Wigner function ( |l4| ) and applying the Liouville 
equation dp/dr — we find: 

dn 8 (x, P ) _^df a ^(x,p) _ f ^ ^ ipy „ d Da , 



0x» =Jdve-^<^ix,y)>. (16) 

The substitution of the equation of motion (^) leads to the relation 

""^/""'(s.P) = -i J dye" 1 ™ < [P^ (x,y), H(r)] > . (17) 

Let us assume that the vector p* is time-like and hence can fix the corresponding direction in 
the Minkowsky space, i.e. 

n M d g U M = p »/^ y u 2 = 1. (1 8 ) 

Substitution of relation (|l8|) into Eq. (^) leads to 

Plt aHr ,P (x,p) = -iy/? I dye- 1 ™ < [P a >\x,y),H{r)] > . (19) 



The right-hand side of this equation is the Vlasov-like drift integral. This general result for the 
mean field approximation holds for different internal degrees of freedom (spin, isospin, color, 
flavor ...). In the following we will drop the spinor indices a,/3 in the notation. 



3. Application of the method to the Walecka model 
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3.1. General structure of the transport equation 

The simplest version of quantum hadro dynamics is the Walecka model of relativistic nuclear 
matter E3. The Lagrange density for the nucleon (tp), the neutral scalar ((f)) and the vector 
(oj m ) mesonic fields reads: 

= \i> 3 ) ,/^M#+^( 5s fj„^ 7 >+i[a^cf^my]if''''^+^^ 2 , (20) 

where M, m s , m v are the masses of the nucleon, scalar and vector mesons, respectively; g s , g v 
are the coupling constants, F^ v = d^u) v — d"^ '. The meson fields are approximated by their 
mean fields, given by <f> =< <f> > and lo^ —< w M > where the symbol < ... >= Tr{...p} denotes 
the procedure of statistical averaging with the density matrix p in the Heisenberg picture. 

A concrete definition of the drift integral is feasible only by specifying the Hamiltonian of 
the system (Q). In the considered model, the Hamiltonian ( |2l[ ) is a bilinear combination of 
the field operators ip(x) and 4>(x), consequently the truncation problem does not arise here. 
For the Walecka model ( pc) ) with the energy- momentum tensor we obtain 

H(t) = -J da(n)^C-Y 8^ -AT - ff.'AvJtf, (21) 
lP(r) = l -J da(n)n A ^{ 7 " d x + 7 A d"}^, (22) 

where 

M* = M - g s c/> 

is the effective nucleon mass in mean field approximation. 

The intermediate calculations of the drift integral in the Walecka model are shown in 
Appendix A. As it ought to be, the resulting relativistic kinetic equation of the Vlasov type 
turns out to be of non-lecal (and non-Markovian) type. The restriction to minimal orders of 
the gradient expansionsEj allows us to write the relativistic Vlasov equation in a local form, 

P^f(x, P) + i^A/*{f, d£f(x, P)} + g v P»F^d" P f(x, P)+ 

\\W,dlf{x,P)]+iM*[rj{x,P)]-±g v F^[W,d» P f{x,P)] = 0, (23) 

whspeiP = p—g v u is the kinetic momentum. The derivation of this kinetic equation is reviewed 
inE3o. Although the model ( po[ ) is not gauge-invariant, this property emerges when m v —* 0, 
therefore in the Walsclta model it is convenient to use the gauge-invariant generalization of the 
Wigner functiorMEaEl In the case of small mesonic field gradients th is g eneralization can be 
done by the replacement p — > P — p — g v tu in (^), this is fulfilled in (p3|). 

The suggested method of the derivation of the relativistic kinetic equation of the Vlasov 
type does not lead to the mass shell condition itself. However the space-like Eqs. of motion 
(g) generate also a set of equations for the Wigner function. Indeed, according to Ech-i(|^) the 
result of the operation (x) ^ on the Wigner function (jlij) leads to the equationsu 

■£zf{x, P ^- % J d y e ~ ,pv < [ p ( x = »). nM ( r )i > • ( 24 ) 

The absence of an additional restriction to the Wigner function is not obvious in general. 
However, it is easy to show that these equations do not contain new information in the mean 
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field approximation of the Walecka model. The operator II M (r) is defined by relation ( p2| ) 
which does not contain the mass and the mean field dependencies and, consequently, the 
Eq. ( p3| ) can not result in the mass shell condition. This conclusion is confirmed by direct 
calculations of the right part of Eq. (^) with the operator (psj) 

afr/fo p ) + ^MtV, u x dif{x, p)] + i^X^fix, p) - Mt'V, f{x, p)}) = o. (25) 

One could assume that the mass shell condition is contained inside of the kinetic equation 
However it will be shown in the following section that this assumption can't be confirmed. 

3.2. Dirac decomposition and properties of the relativistic Vlasov equation 

In order to verify the physical consistency of the relativistic kinetic equation of the Vlasov 
type (^), let us consider the simplest case of the spin saturated system where the spin- 
dependent effects can be neglected. Since the pseudo-scalar and the pseudo-vector contribu- 
tions in the modsl-|(bp|) are absent, the Wigner function decomposition on the basis of the 
Clifford algebra isO'ET 

f(x, P) = / S (x, P) + /J (x, pyf + f%, (x, P)o*» , (26) 

where o» v = i{ 7 ", 7 "} an d 

/ S = iTr{/}, fl = iTr{y/}, /£ = ^Tr{^ f} , (27) 

The symbol Tr denotes the trace with respect to the spinor indices. Substitution of the relation 
( p6[ ) into the kinetic equation ( p^ ) leads to the system of equations for the decomposition 
coefficients (§?])■ Let us write this system in the quasi-classical limit Ti — > when /J„ = 0. 
Projection of the Dirac- scalar and vector part of the Wigner function leads to 



P^f S (x,P)+d x v M*P»d P ft(x,P)+g v P»F^d P f S {x,P) = (28) 
tfl(x,P) +P x dZrf(x,P) - P»dZfi(x,P) +dZM*P x d P f S {x,P) + 
g v P^F^d P fY(x,P)-g v P x F^d^f^(x,P)+g v P^F x ,d P f^(x,P) = (29) 



and 



P ll tf(x,P) = P u fi{x,P). (30) 

For the description of a spin saturated system it is appropriate to find the closed relativistic 
kinetic equation for the scalar part of the Wigner function f s (x,P). 

After a few algebraic transformations, we end up with the following kinetic equation for 
the scalar Wigner function: 

iS 



P l tdZf(x,P)+VP*(dZM*)dp b (x,P)+g v pi*F li „d P f(x,P)=0. (31) 

Let us now discuss some feature of this relativistic Vlasov equation. Firstly, the kinetic 
equation (J^) occurs without the mass shell condition: this condition can be introduced only 
from outside as an additional restriction. On the mass-shell 

P 2 = M* 2 , (32) 
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Eq. (|3lj) turns into the |^rtrk-jW vn relativistic kinetic equation of the Vlasov type for spin 
saturated nuclear matteraoaB'EZI 

P fl dZf s (x,P)+M*d%M*dZf s (x,P)+g v P>*F flv d v P f s (x,P)=0. (33) 

The derived kinetic equation fulfills the conversation laws. For example it is easy to check 
with help of ( |3l| ) that the baryon current density 

f(x) = / dPTc{^f(x,P)} = 4 / dPf£(x,P) = 4 fdpL^lp^ (34) 



P 2 

is fulfilled by the continuity equation 

divj(aO=0. (35) 

The entropy conservation is verified analogously and holds in the investigated mean field 
approximation (no collisions). 

4. Application of the method to the Nambu Jona-Lasinio-model 

In this section we want to derive a relativistic kinetic equation of the Vlasov type for NJL 
model ELiThe NJL model was successfully applied for the description of chiral symmetry 
breaking E3. The effective Lagrangian reads 



£(x) = ^(x)<9 M 7 M ?/> - m,( l ^(x)ip{x) - G 



(ip(x)tp(x)) 2 + {^{x)i^ 5 ij}{x)y 



(36) 



where mo is the current quark mass and G the coupling constant. Within this model study 
we neglect color and flavor degrees of freedom. We restrict to the mean field approximation 
(Hartree approximation) and obtain 

^Hartrcc^-j _ ^^d^^ tp(x) — mo?/> (x)ifj (x) — a (x)ip( x)%j) (x) — n (x)4>(x)i'y 5 (x) , (37) 

where the mean fields in the scalar and pseudoscalar channel are given as 

a(x) = -G < i>(x)ip(x) > , (38) 
tt(x) = -G < tp{x)i-y 5 ^(x) > . (39) 

The corresponding Hamilton density reads 

H int (x) = i>(x)4>{x)ip{x) , (40) 

with the channel decomposition 

4>al3(x) = CT{x)5 a p + n{x)i^ S af3 . (41) 

The substitution of Eq. ( E o j) into the kinetic equation (^9|) leads to the following equation for 
the model Lagrangian (13 71) in lowest order of the gradient expansion 



5> 

1 



P^f(x,p) + -[frj* 1 ,dlf(x,p)] + im Q \p~, f(x,p)] 



+ iW>(x)f(x,p) - f{x,p)4>{x)p} + ~{pXd x a <f>)d£f(x,p) + (d»f(x,p)(d^)p} = 0. (42) 
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The kinetic part of this equation equals to the corresponding terms in Eq. (£3|) . The last two 
terms are result of the interaction within the NJL model. 

In order to derive a closed equation for a spin saturated system we use the Clifford decom- 
position of the Wigner function 

/ = r+ns.r+7 M /;+7' i 7 5 /M A . (43) 

In the quasiparticle approximation we can extract the following set of equations for the coef- 
ficients of the decomposition (Eq. (^)) 

P&f + (dla)p"d$fZ 
Ppd%f p - (d*n)p"d£ff 
P^fx +p\d£f% -p^dlfl +p\{{d^o-)d%f a - {d*n)d£ f} 

M a f + nf 
PuL -PvJv 

where M a = mo + o- In order to proceed we will restrict to the chiral limit (mo = 0) and 
introduce the chiral invariant mass A4 = \/ a 2 + n 2 and the transformed Wigner function as 
T = f s /a. On the mass shell (p 2 = M 2 ) we obtain for the scalar part (JF) 

p„d»F(x,p) + M(x)d*M(x)d£F(x,p) = (49) 

using the additional relation pd x o~ = pd x ir = 0. 

Thei relativistic equation of the Vlasov type (^9|) for the NJL model has the same form 
as in Ej. Starting from this equation it is very interesting to study the dependence of the 
scalar Wigner function (quark condensate) as a function of space and timeE3. Therefore it is 
obviously necessary to define reasonable initial conditions. Such a numerical investigation is 
in progress but out of the scope of this work. 






(44) 





(45) 





(46) 





(47) 


o , 


(48) 


(mo 


= 0) and 



5. Summary 

We have presented a derivation of a relativistic quantum kinetic equation for a system with 
many internal degrees of freedom. For illustrative simplicity, in our analysis we employed a 
mean field approximation. Our approach is based on a non-equilibrium statistical operator 
and can be considered as a relativistic generalization of the Zubarev method. We obtain a 
transport equation which fulfills the conservation laws. 

We applied the approach to the Walecka model and the NJL model. In particular we 
performed the spinor decomposition of the transport equation for both cases. 

In future studies the Vlasov equation should be improved by the inclusion of collision 
integrals. The solution of the relativistic Vlasov equation for the NJL model is of particular 
interest for thereiuark hadron phase transition, where a quark condensate should evolve with 
time and space£3. An interesting application of the approach developed herein would be the 
numerical study of the Vlasov equation in a microscopic quark model. 
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Appendix A: Details of the derivation of the relativistic kinetic equation of 
the Vlasov type 

In this appendix we add some details of the derivation of Eq. (]23|). The starting equation 
is the generalized kinetic equation (M). It is convenient to write this equation in the form 



pd(x)f(x,p) = -I(x,p), (A.l) 
where the drift integral was introduced 

I(x,p) = iy/ifi [dye- 1 ™ < [P(x,y),H(r)] > . (A.2) 



In the considered model (^), the Hamilton operator is defined by formula (p2]). After calcu- 
lation of the commutator in the relation (A.2) by the rules ([n]) and (|l4|), we get 

/ = / Cm) +/ W +J (/) ) ( A _ 3 ) 
where 7 (m) is the mass part of the drift integral, 

I { ^\x,p) = iM\P,f{x,p)] <xB , (A.4) 

I^ k ' is the kinetic part, 

l{ c!$( x >P) = tfs J dye~ lpv < ^i3(x + )d^ {x„)^ s (x-)p' al +d^ {x + )^„X x +)i><*{ x -)fal3 > > ( A -5) 
and I^*' is the mean field part, 

lj${x,p) = i fdye-^x 



x {i> al Y lS (x-) < tl) l3 (x+)'tps(x-) - i>sj3T 1 s{x+) < ip^(x+)tp a (x-) >},(A.6) 

and 

r«p(x) = -g s <t>{x)8 al3 + gyuj^x)^ . (A.7) 

We use the short notation, x± = x ± y/2. The mass part (A.4) has been written already with 
Wigner functions. In order to express the remaining parts of the drift integral (A. 3) in terms 
of the Wigner function it is necessary to apply the transformation to formulae (A. 5) and (A. 6). 
In formula (A. 5) let us consider the first term from the right part 

dye~ lpy < rpp(x+)d^ (x-)ips(x-) > i> al 

- \ dye~* pv {< $p(x+)d£(x)$s{x-) + ii)i3[x+)d£(x-)il)s(x-) >}i> al 



2 
I 
2 
1 

~~ 2 

Here the following identity was taking into account 



dye wy {< ip^x^dp (x)ips{x-) + (<9 M {x)il>p(x+))ip s (x-) >}f> a ^ 
dv(x)fs0(x,p)ft aj - (A.8) 
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according to definition A M „ = g M „ — Py,Pv{p 2 ) 1 ■ The second term of (A. 5) is transformed by 
analogy 

J dye~' py < {d^{x + )i) 1 {x + ))il) OL {x-) > p sf} = ^d^{x)f al (x,p)ps0. (A.9) 

As a result of substitution of (A. 8), (A. 9) into (A. 5) we get one contribution from the relativistic 
Vlasov equation ( |23| ) 

= \d^)[W,f{x,p)U. (A.10) 

The formulae (A. 4) and (A. 10) have the local form already. The mean field part (A. 6) can be 
written in non-local form relatively to the Wigner function without effort. However, in order 
to write this fragment of the drift integral in local form, we perform a gradient expansionuj) 
of the meson fields (A. 7) with respect to the rapid variable 

T a p[x±) ~ Y a p{x) ± ^y ll d' 1 {x)Y af3 (x) . 
The substitution of the decomposition into (A. 6) leads to the result 

I$j(x,p) - {pT(x)f(x,p) - f{x,p)T(x)i>} aP + ±{^(x)T(x)dMf(x,p) 

+ d ll {p)f{x lV )d^{x)Y{x)p} a 0. (A.ll) 

Taking into account the definition (A. 7), fulfilling the necessary commutations and substituting 
the canonical momentum p^ by the kinetic one P^ — p M — g v u)^ , the relation (A.ll) can be 
rewritten in the form of the Eq. (p^). 

References 

[I] QUARK MATTER '96, Nucl. Phys. A 610 (1996). 

[2] S.R. de Groot, WA. van Leeuwen and Ch. G. van Weert, Relativistic Kinetic Theory (North 

- Holland Publishing Company, 1980). 
[3] H.W. Barz, B. Kampfer and B. Lukacs, Sov. J. Part. Nuclei 18, 1234 (1987). 
[4] G. Ropke et al., Sov. J. Part. Nuclei 21, 364 (1990). 
[5] V.D. Toneev et al., Sov. J. Part. Nuclei 17, 1093 (1986). 
[6] B. Blattel, V. Koch and U. Mosel, Rep. Prog. Phys. 56, 1 (1993). 
[7] W. Botermans and R. Malfliet, Phys. Rep. 198, 115 (1990). 
[8] W. Cassing and U. Mosel, Prog. Part. Nucl. Phys. 25, 235 (1990). 
[9] PA. Henning, Phys. Rep. 253, 1 (1995). 

[10] D.N. Voskresensky et al., Int. J. of Modern Phys. E 4, 1 (1995). 

[II] H.-Th. Elze and U. Heinz, Phys. Rep. 183, 81 (1989). 

[12] L.P. Kadanoff and G.Baym, Quantum Statistical Mechanics (Benjamin, 1962). 

[13] PC. Martin and J. Schwinger, Phys. Rev. 115, 1342 (1959). 

[14] L.V. Keldysh, Sov. Phys. JETP 20, 1018 (1964). 

[15] S. Mrowczyriski and U. Heinz, Ann. Phys. (USA) 229, 1 (1994). 

[16] Y. Kluger et al., Phys. Rev. D 45, 4659 (1992). 

[17] D.N. Zubarev, Nonequilibrium Statistical Thermodynamics (Plenum Press, 1974). 
[18] D.N. Zubarev, V.G. Morozov and G. Ropke Statistical Mechanics of Nonequilibrium Pro- 
cesses (Akademie Verlag, 1996). 
[19] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 



Relativistic quantum kinetic equation. 



[20] W.-M. Zhang and L. Wilets, Phys. Rev. C 45, 1900 (1992). 

[21] V.G. Kadyshevsky, Sov. Phys.-JETP 46, 654 (1964); J. Theor. Math. Phys. 32, 34 (1977). 
[22] B.L. Altschuler and A.O. Bravinsky, Uspekhi Fisicheskikh Nauk (Russia) 166, 459 (1996). 
[23] G. Maino ct al., Nuovo Cim. 109 A, 1699 (1996). 

[24] S.A. Smolyansky et al., Proc. XIII Int. Sem. on High Energy Physics Problems, Dubna, 

Russia, 2-7 Sept., 1996, in press. 
[25] H.-Th. Elze, M. Gyulassy and D. Vasak, Nucl. Phys. B 276, 706 (1986). 
[26] S. Fujita, Introduction to Non - Equilibrium Statistical Mechanics (Saunders, 1966). 
[27] H.-Th. Elze et al., Mod. Phys. Lett. A 2, 451 (1987). 
[28] P.T. Mathews, Phys. Rev. 75, 1270 (1949); N.M. Kroll, ibid, 1321. 
[29] Y. Nambu, G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 
[30] U. Vogl, W. Weise, Progr. Part. Nucl. Phys. 27, 195 (1991); 

M.K. Volkov, Ann. Phys. (NY) 159, 282 (1984); 

S.P. Klevansky, Rev. Mod. Phys. 64, 649 (1992); 

T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 
[31] W. Florkowski et al., Ann. Phys. (NY), 245, 445 (1996). 
[32] I.N. Mishustin and O. Scavenius, Phys. Lett. B 396, 33 (1997). 



